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Radiation from relativistic jets in blazars and the efficient 
dissipation of their bulk energy via photon breeding 
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ABSTRACT 

High-energy photons propagating in the magnetised medium with large velocity gradients can 
mediate energy and momentum exchange. Conversion of these photons into electron-positron 
pairs in the field of soft photons with the consequent isotropization and emission of new high- 
energy photons by Compton scattering can lead to the runaway cascade of the high-energy 
photons and electron-positron pairs fed by the bulk energy of the flow. This is the essence of 
the photon breeding mechanism. 

We study the problem of high-energy emission of relativistic jets in blazars via photon 
breeding mechanism using 2D ballistic model for the jet with the detailed treatment of particle 
propagation and interactions. Our numerical simulations from first principles demonstrate 
that a jet propagating in the soft radiation field of broad emission-line region can convert a 
significant fraction (up to 80 per cent) of its total power into radiation. 

We show that the gamma-ray background of similar energy density as observed at Earth 
is sufficient to trigger the photon breeding. The considered mechanism produces a population 
of high-energy leptons and, therefore, alleviates the need for Fermi-type particle acceleration 
models in relativistic flows. The mechanism reproduces basic spectral features observed in 
blazars including the blazar sequence (shift of spectral peaks towards lower energies with 
increasing luminosity). The significant deceleration of the jet at subparsec scales and the 
transversal gradient of the Lorentz factor (so called structured jet) predicted by the model 
reconcile the discrepancy between the high Doppler factors determined by the fits to the spec- 
tra of TeV blazars and the low apparent velocities observed at VLBI scales. The mechanism 
produces significantly broader angular distribution of radiation than that predicted by a simple 
model assuming the isotropic emission in the jet frame. This helps to reconcile the observed 
statistics and luminosity ratio of FR I and BL Lac objects with the large Lorentz factors of the 
jets as well as to explain the high level of the TeV emission in the radio galaxy M87. 

We also discuss other possible sites of operation of the photon breeding mechanism and 
demonstrate that the accretion disc radiation at the scale of about 100 Schwarzschild radii and 
the infrared radiation from the dust at a parsec scale can serve as targets for photon breeding. 

Key words: acceleration of particles - galaxies: active - BL Lacertae objects: general - 
galaxies: jets - gamma-rays: theory - radiation mechanisms: nonthermal 



1 INTRODUCTION 

The large ga mma-ray isotropic lumi nosities of the blazars 10 47 - 
10 49 erg s _1 dMukheriee et al.|[l997t) require a very powerful en- 
ergy source. There is a general consensus that this emission is pro- 
duced by relativistic jets emanating from the vicinity of a black 
hole. The beaming of radiation into a solid angle Q ~ 0.01 sr re- 
duces the total luminosity to 10 45 -10 46 erg s _1 , which is still com- 
parable to the bolometric luminosity of a luminous quasar. This im- 
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plies a high efficiency of the bulk energy conversion into gamma- 
ray em ission, because the accretion luminosity is similar to the jet 
power dRawlings & Saunderslll99lh . 

What kind of energy is available for the high-energy radia- 
tion of blazars: the jet internal energy or the total energy carried 
by the jet through the external environmen t? The first scenario can 
be as sociated with the internal shocks (e.gjRees 1978; Spad a et ail 
l200ll) or magnetic reconnection (see e.g. lLvutikovll2003l) . Then the 
efficiency is limited by the dispersion of the jet Lorentz factor or 
by the share of magnetic field which can reconnect, which depends 
on the field structure. 
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The second scenario provides a larger potential energy budget, 
but requires a "friction" between the jet and the external environ- 
ment. A promising mechanism for such a friction can be associ- 
ated with the exchange of neutral particles between the two me- 
dia which move with respect to each other. The simplest version 
of s uch scenario is the exch ange of photons by Compton scatter- 
ing dArav & Bege lman 1992). If the energy flux in the quasar jet is 
dominated by cold protons then the Thomson optical depth across 
the jet (assuming one electron per proton) is 



MR) 



Li 



Rj 27rm p c 3 rj 



2.3 x 10" 



RneY- s ' 



(1) 



where Lj is the jet power, _Rj is the radius of the jet cross-section, 
R is the distance from the black hole, 9 = RJ R and Tj are the 
jet opening angle and the initial Lorentz factory Such a small opti- 
cal depth is insufficient to provide a considerable photon viscosity. 
However, the optical depth for high-energy photons against pair 
production on the low-energy photons can be high enough. For ex- 
ample, a luminous accretion disc of a quasar emitting soft radia- 
tion with luminosity Ld ~ 10 45 erg s _1 at a typical temperature 
8 = kT/m e c 2 = 10~ 5 , provides a high opacity across the jet (see 
IStern & Po utanen 2006, hereafter SP06): 



t t1 (R) = 60 J*g-(1W), 



(2) 



whic h remains more than u nity up to the pars ec scale. 

iDerishev et alj ( 120030 and IStera d2003h suggested a mecha- 
nism of dissipation of the bulk energy of relativistic fluids through 
the particle exchange with charge conversion (the converter mech- 
anism). When this mechanism supplies energy to e pairs and pho- 
tons, it can work in a runaway manner similar to a chain reaction 
in a supercritical nuclear pile. Here the role of breeding neutrons 
belongs to the high-energy photons, which extract energy from the 
fluid and breed exponentially. In this way the converter mechanism 
takes the fo rm of a sup ercritical photon breeding (hereafter, photon 
breeding). ISterrJ J2003I) studied numerically the operation of the 
photon breeding mechanism in ultrarelativistic shocks propagating 
in a moderately dense environment, showing the high efficiency of 
the bulk energy dissipation. A similar problem for relativistic jets in 
active galactic nuclei (AGNs) was investigated by us (SP06), using 
a simple ID ballistic model of the jet and a detailed treatment of 
particle propagation and interactions. We showed that up to 20 per 
cent of the jet kinetic energy can be converted into gamma-rays, 
if the following conditions are met: Lorentz factor above ~ 4, a 
sharp transition layer between the jet and the external environment, 
the presence of the ambient soft photon field, and a transversal or a 
chaotic, not very strong, magnetic field. 

The present paper is a development of the approach presented 
in SP06. The main improvement is 2D treatment of the fluid (in- 
stead of ID in SP06) and a better numerical resolution. The 2D 
treatment allows us to track the evolution of the system for much 
longer time, when it reaches a steady state. We also study the jet 
emission at a wider range of distances and with different sources of 
the soft photon background. We still neglect the internal pressure 
and treat the jet dynamics in a ballistic approximation. 

In Section[2]we describe our model of the jet, the external soft 
radiation and the scheme of the numerical simulation. In Section 
[3] we describe the parameters of our simulations and present the 



1 We use standard notations Q = 10 X Q X in cgs units and for dimension- 
less variables. 



results: the light curves of escaping photons, their spectra, distri- 
butions of the Lorentz factor of the decelerating jet, and the energy 
distribution of pairs produced in the jet. The properties of the model 
and its limitations are discussed in Section[4] In Section|5]we dis- 
cuss the astrophysical implications of our model. We conclude in 
Section|6] 



2 PROBLEM SETUP 

2.1 Photon breeding mechanism 

Let us first summarise the basic physics of the photon breeding 
mechanism in the context of relativistic jets (see SP06 for more 
details). Consider a jet with bulk Lorentz factor Tj propagating 
through the stationary medium. Assume that the medium is filled 
with the background soft photons [e.g. from the accretion disc, 
broad-emission line region (BLR), dusty torus, etc.] and some seed 
high-energy photons (e.g. extragalactic gamma-ray background). 
An external high-energy photon of energy e enters the jet and inter- 
acts with a soft (disc) photon producing an electron-positron pair0 
The typical energy of the high-energy photon interacting with the 
blackbody radiation field of temperature O is a few times Q~ . The 
time-averaged Lorentz factor of the produced pair (as measured in 
the external frame) gyrating in the magnetic field of the jet, be- 
comes 7 ~ rfe ~ r|G -1 . The electrons and positrons in the jet 
Comptonize soft photons (internal synchrotron or external) up to 
high energies|f] Some of these photons leave the jet and produce 
pairs in the external environment, which gyrate in the magnetic 
field and Comptonize soft photons more or less isotropically. Some 
of these Comptonized high-energy photons enter the jet again. 

In this cycle, the energy gain, ~T?, is provided by the 
isotropization of the charged particles in the jet frame and is taken 
from the bulk energy of the flow. Other steps in the cycle are the 
energy sinks. The whole process proceeds in a runaway regime, 
with the total energy in photons and relativistic particles increasing 
exponentially, if the amplification coefficient (energy gain in one 
cycle) is larger than unity. As shown in SP06, the amplification co- 
efficient is A < r?/10, implying the minimum Lorentz factor of 
the flow needed for the process to operate Tj ~4. This theoretical 
minimum is, however, difficult to reach as other conditions should 
be met. 

The preferred conditions are: the weak jet magnetic field 
(which reduces the energy sink by synchrotron radiation) and large 
energy density of the external isotropic radiation in the jet frame 
(which increases the Compton losses). The latter requirement trans- 
lates into high Tj and/or high fraction of the disc radiation scat- 
tered/reprocessed in the ambient medium. A sufficiently sharp 
boundary between the jet and the external environment also in- 
creases the efficiency of the process; however, as we show in this 
paper, it is not a necessary condition. 



2 The photon energies, in units of the electron rest mass m e c 2 , are denoted 
as x and e for photons of low (< 1) and high ( > 1) energies, respectively. 
We define also the power-law index a = — d log F(x)/d log x, where 
F(x) is the photon energy flux. 

3 Depending on the parameters the seed soft photons are either from the 
external medium or the synchrotron photons produced in the jet. In the first 
case, the inverse Compton (IC) scattered radiation is called the external radi- 
ation Compton (ERC) and in the second case, the synchrotron self-Compton 
(SSC). 
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2.2 Model of the jet 



2.3.2 External radiation 



We consider a jet propagating from the centre of the quasar accre- 
tion disc of luminosity Ld through the medium filled with the soft 
photons from the disc as well as scattered in the environment. The 
total jet power scales with the disc luminosity Lj — tjjLd. The 
comoving value of the magnetic field in the jet Bj (its direction is 
transversal by assumption) is related to the Poynting flux carried by 
the (two-sided) jet 

B 2 

Lb = mLj = ■^ r 27r J R J 2 r j 2 c ~ 8xlO 43 B 2 i? 2 7 (0rj) 2 erg s~\(3) 

which is assumed to be independent of distance (i.e. the magnetic 
field scales as Bj oc 1/R). The ratio of the Poynting flux to the 
total jet power ijb affects the ratio of synchrotron and Compton 
losses. The form of the rest of the energy content is not important 
in our approach as we ignore the internal pressure. We just assume 
that the jet does not contain high-energy interacting particles except 
those that are produced by the studied mechanism. If the rest of 
the jet energy is dominated by cold protons, the jet kinetic power 
(including the rest mass) is 



L k = (1 - r]B)Lj = Fj Mc = 27r(i?6>) 2 rjm p c 3 n p , 



(4) 



where M is the mass outflow rate and n p is the proton concentra- 
tion (in lab frame). The total jet power is then Lj = Lk + Lb. 



2.3 External environment 

2.3.1 Accretion disc radiation 

As discussed in SP06, the photon breeding mechanism requires 
the existence of the external soft radiation field which is opaque 
for high-energy photons against pair production. The emission of 
the accretion disc is opaque for high-energy (~ 10 5 MeV) pho- 
tons moving across the jet up to the distance of a few parsec for 
the disc luminosity Ld ~ 10 45 erg s _1 . We characterise this ra- 
diation by the multicolour disc spectrum assuming a power-law 
T{w) = Tmaxfro/roin) -3 / 4 dependence of temperature on the 
disc radius w, with the ratio of the outer to inner disc radius 
a^out/win = 10 4 . Here T max = 5 eV is the maximum disc tem- 
perature at OTi n = 3i?s and Rs = 2GM/c 2 is the black hole 
Schwarzschild radius. The dimensionless maximal temperature is 
max = fcT max /m e c 2 = 10 -5 . We neglect the expected depen- 
dence on the luminosity T max oc L^ 4 as it is a small effect. 

At distances much larger than the disc size, the disc radiation 
is directed along the jet and does not interact with the gamma-rays 
produced in the jet. Thus the conversion of these gamma-rays into 
electron-positron pairs requires a source of transversal soft photons. 
Below we consider three possibilities for such sources at different 
distance scales. 

A typical AGN in addition to the thermal UV radiation emits 
nonthermal X-ray continuum extending up to ~ 100 keV constitut- 
ing about 10 per cent of total AGN luminosity. This radiation was 
implemented in the model by SP06, however, in the broad emis- 
sion line region (BLR) the X-rays provide very little opacity for 
pair production, and therefore we do not include them in this work 
in order to reduce the number of parameters. Nevertheless, the X- 
rays can be important at smaller distances from the black hole and 
their role is a matter of future studies. 



• Model A: This model represents the isotropic radiation from 
the BLR at a subparsec distance R — 2 x 10 17 cm (as in SP06). 
The soft external radiation in this case is the disc radiation repro- 
cessed and scattered by clouds in the BLR with an admixture of a 
softer infrared (IR) radiation from larger distances. The resulting 
spectrum has a complicated form, which we parametrize here (as 
in SP06) by a cutoff power-law Fblr{x) oc x~ a exp(— a;/a; max ) 
extending from the far-infrared x mln ~ 10~ 9 to the UV band, as 
we assume a; max = max . The normalization of this spectrum is 
given by the parameter rji = Ui so /Ud, which defines the ratio of 
the energy density of the isotropic component to that of the direct 
disc radiation. We use rji = 0.05 in most of the simulations. 

The results are sensitive to the spectral slope a which describes 
the relative number of soft photons. Here we set a = 0.4 which 
implies that the ratio of the isotropic energy density in the IR range 
to that of the direct disc radiation is ~ 0.003 (for rji — 0.05). 
Such IR contribution can be supplied by the dust at a parsec scale 
which can absorb and reemit in the IR range a substantial part of 
the quasar luminosity. 

• Model B: It describes the radiation field at a few parsec scale, 
where the isotropic component is dominated by the IR radiation of 
the dusty torus heated by the disc radiation. The isotropic compo- 
nent of the radiation field is represented as a blackbody with tem- 
perature T = 600 K (6 = 10~ 7 ) and rji = 0.1 or 0.3. 

• Model C: No separate isotropic component is present. 
The radiation field is pro vided by the accretion disc only 
jDermer & Schlickeiserl 19931) . At R ~ 10 16 cm the angular size of 
the accretion disc is large enough to illuminate the jet under large 
angles. We take a standard multicolour disc, but for simplicity as- 
sume that there is one-to-one correspondence between the emitted 
photon energy x and the angle <f> between the jet direction and the 
photon momentum: 



tan <p — — = — — , 

v r r \e„ 



2.3.3 External matter 
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The external environment is kept at rest in spite of the fact that it 
gets a considerable momentum from photons producing pairs out- 
side the jet. This momentum will cause an entrainment of the ex- 
ternal medium which is not easy to simulate. It is important that 
this motion should be nonrelativistic, otherwise the gradient of the 
Lorentz factor would be significantly lower than that in the sim- 
ulations. We discuss the quantitative condition for this in Section 

m 

The magnetic field in the external medium L? cxt is assumed to 
be purely transversal with its direction randomly distributed in the 
plane normal to the jet. The Larmor radius of the highest energy 
pairs, J?L,ie = 0.01578/B cx t.-3, should not exceed the jet radius, 
otherwise the efficiency of the breeding cycle is reduced. Therefore 
we take _B oxt = 10~ 3 G, which is two to three orders of magnitude 
higher than that in the interstellar medium, but still realistic for the 
sub-parsec region of an AGN. 



2.4 Model parameters 

The efficiency of the considered mechanism depends not directly 
on the luminosities, but on the dimensionless compactness param- 
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eters. The compactnesses of the radiation and of the magnetic field 
can be defined as 



t = —^2 U Rj , 



(6) 



where U is the corresponding energy density. The compactness is 
useful when estimating the cooling rate of pairs. The comoving 
magnetic compactness is defined as 



e B = 
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The photon compactness, defined in the external frame, has two 
components: the direct accretion disc radiation 
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(with factor 2 in front of La coming from the Lambert law) and the 
isotropic component £i so — rjild- With these definitions, we can 
write the equation for electron energy losses in the jet frame (for 
model A) as 



dry 



(9) 



where the term with £b describes the synchrotron losses and the 
second and third terms - the Compton losses (in Thomson limit) 
on the isotropic radiation field (with a factor ~ T? coming from 
the Lorentz transformation to the jet frame) and on the direct disc 
radiation, respectively. Because the soft photons can interact with 
high-energy pairs in Klein-Nishina regime, the second term in re- 
ality is much reduced. Equation (|9) does not include SSC losses, 
which are zero at the beginning of the photon breeding cascade (be- 
cause it takes time to built up the synchrotron energy density), but, 
in principle, can be significant at the steady-state regime. The pairs 
produced outside of the jet lose their energy mostly by Compton 
scattering of the disc radiation (the synchrotron losses are negligi- 
ble because of the low magnetic field there) with l& determining 
the cooling rate: 

ft * ~^ 2 - (10) 
The jet compactness 

e ^^^R- C -° A3L ^ (H) 
is related to the Thomson optical depth through the jet 

t p = —-^(1 - Vb) = 2.3 x 1(T 5 (1 - r ?B )L j ,45^ 1 1 6 rr l 1 .(12) 



2.5 Simulation method 

2.5.1 Radiative processes 

The numerical simulation of radiative processes is based on the 
Large Par ticle Monte-Carlo code (LPMC) described in Stern 
Jl985h and lStern et al. The present version of the code han- 

dles Compton scattering, synchrotron radiation, photon-photon pair 
production, and pair annihilation. All these processes are repro- 
duced without any simplifications at the micro-physics level. Syn- 
chrotron self-absorption was neglected as it consumes too much 
computing power and is not important at the considered conditions. 
The number of large particles (LPs) is 2 19 = 524288 for most of 



the simulations (in SP06 we used 2 17 LPs), which is still insuffi- 
cient to completely get rid of numerical fluctuations. A substan- 
tially better resolution will require massive parallel computations 
and a smoothing technique which has to be developed. 



2.5.2 Particle kinematics 

The schem e of particl e tracking in the relativistic fluid was de- 
veloped in ISternl d2003l) and described in details in SP06. The 4- 
momentum of electron/positron LPs are defined in the comoving 
frame of the fluid. If the electron/positron has a critical Lorentz 
factor 
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(13) 



it loses half of its energy in a half of Larmor orbit to synchrotron 
radiation. This requires a fine particle tracking at the first Larmor 
orbit for the comoving Lorentz factor of the particle 7 > 7 max , and, 
therefore, the comoving tracking step is limited to At = 0.17?l/c, 
where 7? L = 1.7 x 10 3 7/Bj cm is the Larmor radius. 

In the external medium, synchrotron losses are small and an 
electron is cooled by the disc radiation. The condition that the elec- 
tron cooling time-scale is longer than the inverse of Larmor fre- 
quency can be written as 



1 2ttR l 

> ~W 

which translates to the limit on the Lorentz factor 

1/2 

. .. f : / ./ 1 i 1 i ; i-V v 1 :: \ 

7 < 7cxt .max 
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(14) 



(15) 



In our simulations, we exactly track the particle motion for the 
duration of the first Larmor orbit for the comoving 7 > 2 x 10 5 . For 
smaller 7 and at further orbits we neglect the dependence between 
the comoving direction and the energy and sample the direction of 
the particle assuming a uniform distribution of its gyration phase in 
the comoving system (see details in SP06). 

The trajectories of electrons and positrons in the magnetic 
field are simulated assuming a random direction of the field in the 
plane normal to the direction of motion in the jet as well as in the 
external matter. The magnetic field correlation length is assumed to 
be larger than the Larmor radius. 



2.5.3 Geometry and jet dynamics 

Our model of the jet is the same as in SP06 except now we use 2D 
representation of the fluid instead of the ID. We consider a piece 
of the jet centred at distance R from the central source and length 
20 Rj, where the jet radius is Rj = R6 and 9 has the meaning of 
the opening angle, while for simplicity we approximate the jet by 
a cylinder of radius Rj . The boundary condition at the inlet is the 
flow of constant Lorentz factor Tj (see Fig.[T|lQ 

We introduce dimensionless coordinates, where the transver- 
sal distance from the jet axis r and the distance from the central 
source z are measured in units of Rj . We use a fixed Eulerian grid 
with fairly rough spacing: 20 cells in the ^-direction and 64 cells 



4 Our model assumes a well-collimated jet with 9 < l/iy If the open- 
ing angle of the jet is much larger, then the efficiency of photon breeding 
drops, because the probability for the high-energy photons to escape the jet 
decreases. 
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z=0 z=20 




rl=1.5 
rl 

Figure 1. Schematic representation of the jet geometry. The jet is repre- 
sented by a cylinder of radius r = 1. The simulation volume of radius 
r = 1.5 includes also the external medium. At the inlet, 2 = 0, the jet 
has Lorentz factor I\ . At the outlet, z = 20, the terminal Lorentz factor Tf 
depends on radius r. 



Interactions between particles change the energy and momen- 
tum of the fluid. The electron-positron pairs born in the process 
become part of the fluid. The momentum and energy exchange be- 
tween photons and each grid cell of the flow is accumulated during 
the time step At — 0.01-Rj / c and the total 4-momentum of the cell 
is updated after that. The momentum is transferred between cells 
along the jet axis. 

The output of the simulations consists of the photons that es- 
cape from the cylindric "simulation volume" with the boundaries 
r — 1.5, z = 0, and z = 20. We also record the energy spectrum 
of photons and pairs in the simulation volume. 



3 RESULTS 



in the r-direction. The step on r is nonuniform: the spacing is finer 
near the jet boundary. The time unit is Rj / c. 

The fluid in each cell i is characterised by its total energy Ei 
and momentum Pi, which define the cell Lorentz factor: 



The flow is treated as one-dimensional along the cylinder axis, be- 
cause the transversal component of P[ never exceeds 10 -2 of the 
total momentum and the transversal displacement of the flow can 
be neglected. 



2.5.4 Statistical representation 

For a good statistical representation of the process, it is necessary 
to have more or less constant density of the LPs over the volume of 
the cylinder. As the energy (and number) density of particles (pairs 
and photons) grows by orders of magnitude from the inlet of the 
cylinder to the outlet, we need to introduce the dependence of sta- 
tistical weights on z. The statistical weight w of a particle of energy 
e is defined through the energy content we — 2™ s, where n is the 
integer part of z/2 for z < 14 and n = 7 at z > 14. The value s 
varies with time in order to keep the total number of LPs approxi- 
mately constant. This weighting scheme is implemented as follows: 
we split a particle into two identical particles with twice less statis- 
tical weight when it crosses one of the planes z = 2, 4, 6... in the 
direction opposite to the jet motion and we kill a particle with prob- 
ability 0.5 and increase its weight by 2, when it crosses the plane 
in the direction of jet motion. The price for a good resolution at the 
region, where the avalanche starts and grows, is a poor resolution at 
the output. For that reason the Poisson noise in output light curves 
and photon spectra is rather high. 



2.5.5 Simulation setup 

At the start of simulations we set constant Lorentz factor Fj over 
the jet volume (r < 1 and < z < 20, see Fig. Q} and assume a 
sharp boundary between the jet and the surrounding medium. We 
also inject a number of high-energy photons with the total energy 
several orders of magnitude smaller than the steady-state luminos- 
ity of the jet. (There is no injection of any additional photons later 
during the simulations and the high-energy population of particles 
is self-supporting.) The external soft-photon radiation field is con- 
stant during the simulation. 



3.1 Simulation runs 

We consider the distance scale of R ~ 10 16 -10 19 cm and the jet 
Lorentz factor Tj ~ 8-40. We vary the luminosity of the accretion 
disc Ld in the interval ~ 10 43 -10 46 erg s _1 and try different ratios 
of the jet to the disc power. First, we present results of three series 
of simulations, where we vary one of the parameters and nine non- 
systematic runs covering a wider area of parameters to demonstrate 
the working area of the mechanism. In the two series we vary the 
disc and jet luminosity, while in the third series we vary the jet 
Lorentz factor. Most of our simulation runs concern model A. Three 
runs correspond to different astrophysical scenarios described by 
models B and C, where the photon breeding operates. 

(i) Runs 01-05: This series represents magnetically dominated 
jet (i.e. r/B = 1) with the jet power equal to the disc luminosity 
(r/j = 1), which is varied between Ld — 3 x 10 45 erg s _1 down to 
3 x 10 43 erg s _1 , where the photon breeding becomes sub-critical 
and the process does not work (at R — 2 x 10 17 ). 

(ii) Runs 10-16: Similar to the previous series, but for the matter 
dominated jet with rye = 0.2 and r/j — 1. We vary Ld in a wider 
range, because the process for smaller t]b has a lower luminosity 
threshold. 

(iii) Runs 21-25: In this series we have fixed the disc luminosity 
to Ld = 3 x 10 44 erg s _1 use r/j — 1 and r/B = 0.2, and vary the 
initial Lorentz factor of the jet. 

(iv) Runs 31-35: Run 30 demonstrates a photon breeding pro- 
cess starting from the extragalactic gamma-ray background and re- 
sulting in 20 orders of magnitude growth of the high-energy photon 
population (see Section |4~5l >. Runs 31 represents a strongly non- 
linear system, where the jet power exceeds the disc luminosity. Run 
32 is the case of the weakest AGN (Ld = 2 x 10 43 erg s" 1 ) tried in 
this work for model A. Run 33 demonstrates that a high efficiency 
can be achieved at Tj = 10 if the magnetic field is weak. Run 34 is 
a trial with the lowest Lorentz factor, Tj — 8, and small Poynting 
flux. Run 35 is similar to run 12, but for the harder spectrum of 
isotropic external photons (a = instead of a = 0.4). 

(v) Runs AI-A2: They correspond to Model B, where we simu- 
late interaction of a strong jet Lj = 2.3 x 10 45 erg s _1 of Lorentz 
factor Tj = 20 with the IR radiation (6 = 10~ 7 ) of the dusty torus 
at a parsec scale. Run 43 describes Model C, where a relatively 
weak, 77j = 0.1, jet interacts with the direct radiation from the ac- 
cretion disc of Ld = 8 x 10 43 erg s _1 at a distance of a hundred 
gravitational radii (for a 10 8 solar mass black hole). 

The parameters for these runs are summarised in Table 1 . 
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Table 1. Model parameters and results. 



Run 


R a 
cm 


I? 


9° 


r d 

erg s — 1 


ER e 




Vi s 






V 2 f- ■? 

1 j *-iso 




B< 
G 


m 

r P 


tpb 


r f ° 


£off p 


01 


2E17 


20 


0.05 


3E45 


A 


0.05 


1.0 


1.0 


3.3E-3 


6.6E-2 


3E-3 


3.2 





0.5 


18.0 


0.25 


02 


2E17 


20 


0.05 


1E45 


A 


0.05 


1.0 


1.0 


1.1E-3 


2.2E-2 


1E-3 


1.8 





1.1 


18.0 


0.32 


03 


2E17 


20 


0.05 


3E44 


A 


0.05 


1.0 


1.0 


3.3E-4 


6.6E-3 


3E-4 


1.1 





3.7 


18.6 


0.26 


04 


2E17 


20 


0.05 


1E44 


A 


0.05 


1.0 


1.0 


1.1E-4 


2.2E-3 


1E-4 


0.6 





6.2 


19.6 


0.05 


05 


2E17 


20 


0.05 


3E43 


A 


0.05 


1.0 


1.0 


3.3E-5 


6.6E-4 


3E-5 


0.35 





~ 


20.0 





10 


2E17 


20 


0.05 


1E46 


A 


0.05 


1.0 


0.2 


1.1E-2 


2.2E-1 


2E-3 


2.6 


9.3E-5 


0.1 


16.2 


0.42 


11 


2E17 


20 


0.05 


3E45 


A 


0.05 


1.0 


0.2 


3.3E-3 


6.6E-2 


6E-4 


1.4 


2.8E-5 


0.2 


16.2 


0.44 


12 


2E17 


20 


0.05 


1E45 


A 


0.05 


1.0 


0.2 


1.1E-3 


2.2E-2 


2E-4 


0.82 


9.3E-6 


0.42 


14.7 


0.50 


13 


2E17 


20 


0.05 


3E44 


A 


0.05 


1.0 


0.2 


3.3E-4 


6.6E-3 


6E-5 


0.45 


2.8E-6 


1.2 


13.6 


0.56 


14 


2E17 


20 


0.05 


1E44 


A 


0.05 


1.0 


0.2 


1.1E-4 


2.2E-3 


2E-5 


0.26 


9.3E-7 


3.0 


15.6 


0.41 


15 


2E17 


20 


0.05 


5E43 


A 


0.05 


1.0 


0.2 


5.5E-5 


1.1E-3 


1E-5 


0.18 


4.6E-7 


5.2 


18.6 


0.12 


16 


2E17 


20 


0.05 


3E43 


A 


0.05 


1.0 


0.2 


3.3E-5 


6.6E-4 


6E-6 


0.14 


2.8E-7 




20.0 





21 


2E17 


40 


0.05 


3E44 


A 


0.05 


1.0 


0.2 


3.3E-4 


2.6E-3 


1.5E-5 


0.22 


1.4E-6 


0.62 


10.0 


0.82 


22 


2E17 


30 


0.05 


3E44 


A 


0.05 


1.0 


0.2 


3.3E-4 


1.5E-3 


2.7E-5 


0.30 


1.8E-6 


0.78 


10.3 


0.77 


23=13 


2E17 


20 


0.05 


3E44 


A 


0.05 


1.0 


0.2 


3.3E-4 


6.6E-3 


6E-5 


0.45 


2.8E-6 


1.2 


13.6 


0.56 


24 


2E17 


14 


0.05 


3E44 


A 


0.05 


1.0 


0.2 


3.3E-4 


3.2E-3 


1.5E-4 


0.71 


3.9E-6 


3.0 


13.6 


0.14 


25 


2E17 


12 


0.05 


3E44 


A 


0.05 


1.0 


0.2 


3.3E-4 


2.4E-3 


1.9E-4 


0.80 


4.6E-6 




12.0 





30 


2E17 


20 


0.05 


1E44 


A 


0.05 


1.0 


0.12 


1E-4 


2E-3 


1.3E-5 


0.2 


1.0E-6 


2.7 


12.7 


0.53 


31 


2E17 


20 


0.05 


1E45 


A 


0.05 


9.0 


0.13 


1E-3 


2E-2 


1.3E-3 


2 


9.0E-6 


0.95 


17.7 


0.30 


32 


2E17 


20 


0.05 


2E43 


A 


0.05 


1.0 


5E-3 


2.3E-5 


4E-4 


1.3E-7 


0.02 


2.3E-8 


7.0 


19.5 


0.03 


33 


2E17 


10 


0.05 


3E44 


A 


0.05 


0.8 


1.3E-2 


3E-4 


1.5E-3 


1.3E-5 


0.2 


5.9E-6 


1.8 


8.7 


0.31 


34 


2E17 


8 


0.05 


3E44 


A 


0.05 


1.0 


7E-3 


3E-4 


1E-3 


1.3E-5 


0.2 


9.3E-6 


2.6 


7.7 


0.13 


35ssl2 


2E17 


20 


0.05 


1E45 


A8 


0.05 


1.0 


0.2 


1.1E-3 


2.2E-2 


2E-4 


0.82 


9.3E-6 


1.8 


16.2 


0.41 


41 


6E18 


20 


0.033 


1.3E45 


B 


0.3 


1.8 


2E-2 


3E-5 


3.7E-3 


2.6E-6 


0.02 


1.3E-6 


5.4 


18.0 


0.21 


42 


6E18 


20 


0.033 


1.3E45 


B 


0.1 


1.8 


6E-3 


3E-5 


1.2E-3 


6.4E-7 


0.01 


1.3E-6 


6.3 


19.1 


0.14 


43 


6E15 


20 


0.067 


8E43 


C 





1.0 


6E-2 


4E-3 





1.3E-5 


1.0 


2.2E-6 


2.3 


19.2 


0.19 



"Distance from the black hole, 'initial Lorentz factor of the jet. c Half-opening angle of the jet. d Total disc luminosity. e Model of the external radiation 
described in Section l2.3.2l ■' Ratio of the isotropic soft photon energy density to that of the disc. 9 Ratio of the total jet power to the disc luminosity, Lj / Lj. 
^Fraction of the jet power carried by the Poynting flux, Lb/ Lj. 'Disk compactness. ^Compactness of the isotropic external radiation in the jet frame. 
fe Magnetic compactness. 'Magnetic field in the jet. m Thomson optical depth through the jet produced by electrons associated with protons. n e-folding 
time of the photon breeding in Rj /c units. "Terminal Lorentz factor at the jet axis. p Average efficiency of the jet energy conversion into radiation in the 
steady-state. 'Modified model A with a harder spectrum of external isotropic photons (a = 0). 



3.2 Time evolution and the steady-state regime 

The energy release curves for set 2 (runs 10-15) are presented in 
Fig. [2] All curves show an exponential rise at the start with the time 
constant depending on the disc luminosity. The fast exponential rise 
turns to a slower growth when approaching the steady state. A very 
fast rise (time constant ~ 0. lRj /c for run 10) occurs in the case of 
a dense external radiation field when the photon free path length is 
much smaller than the jet radius. Then the photon breeding takes 
place only in a narrow layer at the jet boundary. Such fast breeding 
terminates after a short time, when the thin boundary layer deceler- 
ates. The breeding turns to a slower mode, where the softer photons 
with a longer free path dominate. 

Note that the injection of external high-energy photons takes 
place only at the start of the simulation. Then the regime becomes 
self-supporting. The reason for this is a spatial feedback loop. The 
photon avalanche develops from a smaller to a large z, but some 
high-energy photons produced in the external environment move in 
the opposite direction, from larger z, where there are many high 
energy photons, to smaller z, where they can initiate a new photon 
avalanche. Due to this feedback a high efficiency of emission e c ff = 
Lph/Lj can be reached even if the growth of the photon avalanche 
is slow. 



We can see large fluctuations of the energy release at the sat- 
uration regime (at the stage of the exponential rise, the curves are 
smooth). Fluctuations have the Poisson component which is large 
because of the weighting method: LPs at the outlet have statistical 
weight 2 7 larger than at the inlet. There are also non-Poisson fluc- 
tuation which are caused by fluctuations in the spatial feedback: a 
photon with a large statistical weight moves from the end of the 
active jet fragment to its beginning, interacts there and produces a 
photon avalanche. The number of such photons is small, but the 
amplification coefficient (i.e. the number of secondary photons and 
pairs) can exceed 10 4 . Therefore the number of LPs as large as 
~ 10 6 is still insufficient for a good statistical representation. 

However, statistical noise is not the only reason for the large 
fluctuations in the output photon flux. In runs 10 and 31, where the 
compactness is high, and especially in run 31, where the system is 
strongly non-linear, we see a large amplitude variability at a range 
of time-scales (see Fig.[3](- We performed run 31 with the different 
number of LPs: 2 19 and 2 21 , but the character and the amplitude 
of the fluctuations was similar. Thus, in this case we deal with the 
real instability of a non-linear dynamical system (see Section |4~2l 
for discussion). 
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Figure 2. The power converted into photons (including those born in the 
jet, but not yet escaped) versus time for series of runs 10-15 (from top to 
bottom). 




Time [R/c] 



Figure 3. The escaping photon luminosity in units of the jet power versus 
time for runs 10, 11, 21 and 31. Note that the light curves look different 
compared to Fig. [2] as it takes about 20 time units for photons to escape 
from the simulation cylinder (and the ordinate scale is linear here). 



3.3 Radiative efficiency and the "working area" of the 
process 

The resulting efficiency of the conversion of jet power into radiation 
at the steady-state and the time constant of the exponential photon 
breeding are shown in the last columns of Table 1. Once we fixed 
R, Rj , and the ratio rji in all runs for model A, our main remaining 
parameters are the jet Lorentz factor, the disc luminosity and the 
Poynting flux. The jet power is not important unless it exceeds the 
disc luminosity (as in run 31), making the system strongly non- 
linear. The series of runs 01-05 shows that the process can work 
efficiently in the case of a magnetically dominated jet giving e c a ~ 



1/4. The lowest luminosity when it works is ~ 10 44 erg s _1 . If 
the matter dominates over the magnetic energy by a factor of 5 
(runs 10-16) then the radiative efficiency is higher, above 50 per 
cent, and the working luminosity range slightly extends down to 
5 x 10 43 erg s -1 . The supercritical breeding is still possible at even 
lower disc luminosity, 2 x 10 43 erg s _1 (see run 32), but only if 
the magnetic field is very weak. Below this threshold the photon- 
photon opacity across the jet (see equation[2j becomes insufficient. 
(We note here that the working luminosity limits are given for R = 

2 x 10 17 cm, Tj = 20 and ^ = 0.05.) 

Remarkably, the efficiency at a moderate disc luminosity (~ 

3 x 10 44 erg s _1 , run 13) is higher than for a more luminous quasar 
(runs 10, 11 and 31). The reason is that at a high luminosity the 
spine of the jet is shielded from the high-energy external photons 
by the opaque (respectively to pair production) radiation of the ac- 
cretion disc and synchrotron radiation of pairs in the jet. In runs 
10 and 31, only the photons with energy e < 10 3 can penetrate 
deeply into the jet. Nevertheless, such moderate energy photons 
still can interact with the synchrotron photons (mostly produced in 
the decelerated outer jet layers) inside the jet and support the pho- 
ton breeding. 

As one can expect, the dependence on Fj is very strong. The 
efficiency e a s is close to (Tj — Ti)/Tj, where the average Lorentz 
factor at the outlet TV slightly decreases when Tj increases (see 
Section[33}. 

At a lower Tj the process becomes very sensitive to the ratio 
of the Poynting flux to the isotropic disc luminosity. Indeed, at the 
given magnetic flux Lb the comoving magnetic compactness scales 
as £b oc Le/r? (see Eqs.[7jl, while Compton cooling rate in the 
Thomson limit scales as oc £i BO V?, see equation {9}. Thus the ratio 
of Compton to synchrotron losses scales as r/iVfO 2 /r/B. This ratio 
decreases somewhat because of the Klein-Nishina effect, but it is 
still a strong function on Tj. That is the reason why the runaway 
photon breeding does not happen at Fj < 14 if t)b = 0.2 (series 
21-25). If the magnetic field is very weak, then the process can 
work efficiently at Tj = 10 (run 33) and less efficiently, but still in 
a self-supporting regime at Tj = 8 (run 34). 

All results described above are obtained under an assumption 
that the BLR is relatively abundant with isotropic IR photons, so 
that the photon spectrum can be roughly approximated as a cutoff 
powerlaw with a = 0.4 (see Section |2. 3. 2\ . The spectrum how- 
ever can be harder if, for example, there is no sufficient amount 
of dust around the AGN. Therefore, we made a couple of trials 
for model A with a harder BLR spectrum, a — 0. One trials with 
the same parameters as run 13 (Ld = Lj = 3 X 10 44 erg s~ x ) 
resulted in no self-supporting photon breeding. Another trial with 
(Ld = Lj = 10 45 erg s _1 (run 35 in Table 1) demonstrated pho- 
ton breeding, which was slower than in run 12 (with a — 0.4). 
However, the final steady state was almost the same including the 
radiative efficiency (0.41 versus 0.5 in run 12) and the emitted spec- 
trum. Such a striking difference between the two cases differing 
only by factor 3 in luminosity appears because of the non-linear 
character of the mechanism for high-power jets: soft synchrotron 
photons produced by the jet significantly contribute to the photon- 
photon opacity. The problem with soft photon starvation can be 
resolved if some internal jet activity (e.g. internal shocks) supply 
soft synchrotron photons w hich can be then scattered in the BLR 
dGhisellini & Madau|[l996T) . 

In general, the results of our simulations for distance R from 
the black hole and luminosity L, can be scaled to another distance 
R' if one changes at the same time all luminosities as L' = LR' / R. 
In this case, the compactnesses of all radiation components and of 
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Figure 4. The transversal distribution of the bulk Lorentz factor, (a) Dependence on distance z for run 22 (curves from top to bottom correspond to z = 2 
with step 2 till z = 20). (b) Distribution of the terminal Lorentz factor at the outlet (z = 20) for runs 21-24 corresponding to Tj =40-14. 



the magnetic energy do not change (see equation[6}. Such scaling is 
not exact, because the value of magnetic field changes and the syn- 
chrotron spectra change accordingly. The scaling, however, does 
not affect the criticality of the system as long as the ratio between 
synchrotron and Compton losses remains constant. 

At luminosity smaller than about Ld.min ~ 10 43 i?i7erg s _1 
the opacity for pair-production through the jet becomes smaller 
than unity and the super-criticality condition breaks down. Assum- 
ing the distance to the emission region of R ~ 10 3 i?s, we can 
express the luminosity threshold in terms of the Eddington lumi- 
nosity LcLmin ~ 10" 5 L E dd- 

The distance to the BLR may also depend on the luminos- 
ity of the object. The reverberation mapping gives the relation 
■Rblr.it ~ £^44 JPeterson|[l99l iKaspi et alj|2005l) . This then 

1/2 

gives the condition Ld,44 > 0.1L d 44 , which translates to the lower 
limit on the luminosity Ld > 10 42 erg s _1 , when the photon breed- 
ing can operate. 

The conditions when model A operates can be summarised as 
follows: 

• AtTj > 20 and at Ld, 44 ~ £j,44 > 5J?i7 the photon breeding 
mechanism is efficient and very robust: it works independently on 
the share of magnetic energy in the jet and is not sensitive to the 
spectrum of the external radiation. 

• At Tj > 20 and lower jet power Lj.44 < 5i?i7, when the disc 
luminosity is still sufficiently high, Ld,44 > 0.5i?i7, the mecha- 
nism is sensitive to the spectrum of isotropic radiation in the BLR 
and, if there is a deficit of soft photons, the photon breeding is not 
self-supporting. 

• For Tj < 20 the mechanism becomes sensitive to the magnetic 
field in the jet. The strong field inhibits the photon breeding. When 
the magnetic energy is much smaller that the matter kinetic energy, 
the efficient emission is possible down to Tj ~ 8. 

• The threshold of the disc luminosity is Ld,44 ~ O.I-R17. How- 
ever, the self-supporting emission at such low Ld requires weak 
magnetic field and a soft spectrum of the BLR radiation. 

We have also tried models B and C just to show that the pro- 
cess can work in different situations. Radiative efficiencies of 20 




r=0.99 
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Figure 5. Distribution of the fluid Lorentz factor in z direction for runs 
21-24. Solid curves show the Lorentz factor at the jet axis r = 0, and the 
dashed curves are for the jet boundary at r = 0.99. 

per cent are easily reached. We do not perform a systematic study 
of these models because both scenarios depend on various condi- 
tions that require a different study. Model C, for example, should 
be considered together with the mechanism of the jet launching and 
in model B one should account for the history of the jet flow before 
it reaches the parsec-scale distance. 

3.4 Jet deceleration and the energy release 

Jet deceleration in the ID approach of SP06 took place mainly near 
the jet boundary (r > 0.9), while the centre of the jet was not af- 
fected. In the 2D simulations we can track the deceleration of the 
deeper parts of the jet up to its axis. The transversal distributions of 
the fluid Lorentz factor for runs 21-24 are shown in Fig. [4] Distri- 
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butions of F along the jet axis are shown in Fig. [5] It is remarkable 
that the terminal Lorentz factor T{ depends weakly on the initial 
Fj, actually it even decreases when Tj increases. 

The highest gradient of the Lorentz factor is near the jet 
boundary. In run 22, the Lorentz factor near the boundary drops 
from 30 down to 5 at z = 3 (see Fig.[5}. The jet experiences similar 
deceleration in runs 01, 10, 11 and 31, where the radiation density 
is high. At the outlet the gradient of Tf depends weakly on the dis- 
tance from the axis. It reaches some critical value beyond which 
the photon breeding is not effective anymore. In cases of a lower 
compactness (runs 04, 14, 15, 32) the Lorentz factor changes more 
smoothly with z. At the outlet, dr/dr is still large close to the 
boundary, with the jet getting the form of a fast spine and a slow 
sheath. A large gradient of F means that our ballistic approxima- 
tion rapidly breaks down in the cases of a large compactness or a 
large initial Lorentz factor. We discuss the consequences in Section 

EU 

3.5 Radiation spectra 

Figure [6ja) shows the output spectra of photons emitted in the 
direction of the jet (8 < 10°). The spectra (except those 
for runs 10 and 31 with very high compactness) show a two- 
component synchrotron-Compton structure. Runs 01-03 produce 
a synchrotron peak in the range e ~ 10~ 3 -10 _1 , consistent with 
those observed in BL Lacs, however, the theoretical peaks are less 
prominent that the observed ones (detailed comparison with the ob- 
served spectra is presented in Section |5~Tt , partially due to averag- 
ing over angles. The Compton component at high Ld and Lj (runs 
01,02, 10-12, 31) has a cutoff at 10 4 < e < 10 5 due to the photon- 
photon absorption: the external isotropic BLR radiation provides a 
larger than unity optical depth for the high-energy photons. 

The spectra for a smaller compactness (runs 04, 14, 15, 32, 
as well as runs 41, 42, where the compactness is small because of 
large K), have prominent synchrotron peaks at high energies e ~ 
10°-10 3 . The synchrotron and Compton components are clearly 
separated above the maximal energy of synchrotron radiation: 

e S ,max ~ TnLxSj/Scr « 230Tj , (17) 



which does not depend on the magnetic field (here B CT = 4.4 x 
10 13 G is the critical field). The value of 7 max given by equation 
d 1 3 1 ) is defined as the maximal electron energy (in the jet frame) af- 
ter half the Larmor orbit, i.e. when it turns around from the counter- 
jet direction to the jet propagation direction (see SP06). The derived 
e s ,max assumes that the magnetic field is uniform at the scale of 
the electron Larmor radius and is transversal to the jet. It is not 
sensitive to the soft radiation field because an electron with the co- 
moving Lorentz factor 7 > 10 s interacts with the external radia- 
tion in the deep Klein-Nishina regime and the synchrotron losses 
dominate. Thus, the model predicts a new spectral feature which, 
in principle, can be observed by GLAST. Its detection would be a 
clear signature of the photon breeding mechanism (a distribution 
of diffusively accelerated electrons will hardly extend to the en- 
ergy, where the electron loses a large fraction of its energy in one 
Larmor orbit). Moreover, because of independence of e s , max of the 
magnetic field and other parameters, the Lorentz factor of the jet 
can, in principle, be determined from the detailed fits of the shape 
of the spectral cutoff. 

For the considered model parameters, the dominant source of 
seed soft photon for Compton scattering is BLR. At the start of sim- 
ulations, the SSC losses are small, because the synchrotron energy 
density needs to be built up from zero. In the steady-state regime, 
the ratio of the SSC to the synchrotron losses can be expressed 
as £ e ff?7s,KN/?7B, where ?7s,kn is the fraction of soft synchrotron 
radiation (i.e. Comptonized in Thomson regime) in the total jet lu- 
minosity. We can estimate from Fig. ((Ja) that ry s ,KN < 0.3 for 
runs 01-03 (with 773 = 1) and i) 8 ,kn < 0.1 for runs 10-12 (with 
??b = 0.2). For a lower jet and disc luminosity (runs 03-04, 13-15), 
the electrons in the jet have higher energies and the SSC losses are 
suppressed by the Klein-Nishina effect. Therefore the SSC losses 
are always at least a few times smaller than the synchrotron losses 
and usually are much smaller than the ERC losses. However, the 
Compton losses of pairs in the fast spine of the jet on the syn- 
chrotron photons from the slower sheath can be considerable, at 
least in runs 10 and 31 with the highest compactness. These cases 
can be considered as intermediate between SSC and ERC. 

The radiation spectra at large angles relative to the jet axis 
are shown in Fig.[6jb). We again see two-component, synchrotron 
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Figure 7. Intensity of radiation (per unit solid angle) at different angles to 
the jet axis for run 01 . 



and Compton spectra. The total (angle-integrated) luminosity in the 
gamma-rays emitted at large angles to the jet axis is about two order 
of magnitude smaller than that of the beamed emission, which ap- 
proximately corresponds to the gain factor Vf. The radiation spec- 
tra emitted as a function of the angle 6 between the line of sight and 
the jet axis are shown in Fig. [7] We see that the synchrotron peaks 
shift to lower energies at larger angles and the two components are 
better separated. 

3.6 Angular distributions of radiation 

3.6.1 Radiation from isotropic pairs in the steady jet 

Let us first discuss what kind of angular distributions can be ex- 
pected from a steady (narrow) axisymmetric jet in relativistic mo- 
tion. The observed bolometric luminosity is rel ated to the emitted 
one (in the jet frame) b y the relation (see e.g. iLind & Blandfordl 
ll985l;ISikoraetal.lll997l) 



L{6) 



dL 
dfi 



<5 3 dL' 



(18) 



where 5 = 1/Tj(l — /3j cos#) is the Doppler factor, and df2 = 
d(j>d cos 8. Let us define the amplification factor as the ratio of the 
observed and emitted luminosities 



A{8) 



L{8) 
L'/itr' 



(19) 



If the radiation source is isotropic in the jet frame, then dL' /dQ! = 
L'/4tt, and 



Figure 8. Luminosity (per unit solid angle) as a function of angle between 
the light of sight and the jet axis for run 01 in the three energy bands: IR- 
optical (e < 10 -5 ), X-rays (10~ 3 < e < lO -1 ) and GeV gamma-rays 
(e > 10 3 )(histo grams, arbitrary normalization). The theoretical dependen- 
cies 120) expected for the steady relativistic jet emitting isotropically in the 
comoving frame shown by solid (black) curves (cases with T = 2 and 20). 
Dashed curves represent the synchrotron and SSC emission from pairs in 
2D circular motion, given by equation {24). The solid (red) curve describes 
the jet ERC radiation from isotropic pairs (equation 1221 , while the dotted 
(red) curve is for the jet ERC radiation from 2D pairs (equation 1261 . The 
dot-dashed curve represents the emission of isotropic pairs in the external 
medium Comptonizing radiation beamed along the jet (equation 1291 . The 
emission from 2D pairs in the external medium is shown by long dashes (IC 
of beamed emission, see equation 1281 and short-long dashes (synchrotron 
as well as IC of the isotropic photons, see equation|27j. 



radiation IC scattered by relativistic isotropic electrons (see e.g. 
iRvbicki & LightmarJI 19791) 



dL' 2/3, , 2 



(21) 



with 8' being the angle to the direction of jet propagation. For small 
8, we can rewrite 1 + cos#' ~ S/T, and we get a much sharper 
angular dependence of the ERC radiation (see also IPermeJI 1 9951 
and Fig.(8](: 



A{8) 



3 5^_ 
Iff' 



(22) 



Now the integration over all solid angles shows that the total ob- 
served power is 3/2 larger than the emitted one, which results from 
the non-zero net momentum of the emitted radiation in the jet co- 
moving frame J21t . 



A(P) = TT- 



(20) 



Integrating over angles, we can verify that the total emitted and 
observed powers are the same. Equation d20| l can be applied, for 
example, to the synchrotron radiation in a tangled magnetic field or 
to the SSC radiation by isotropic electrons. 

If soft photons are isotropic in the external medium, they 
are strongly beamed in the direction opposite to jet propagation 
in the jet frame. This results in the angular distribution of the 



3.6.2 Radiation of primary pairs 

In the above we assumed explicitly an isotropic electron distribu- 
tion. However, if the jet magnetic field has a transversal geometry 
(as assumed in this paper), the distribution may significantly devi- 
ate from the isotropic one. For example, the first generation of pairs, 
produced in the jet by the external high-energy photons (which are 
beamed in the jet frame similarly to the external soft photons), has a 
very narrow distribution of pitch angles around n/2. Let us now de- 
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rive the angular distribution of radiation expected from these "2D" 
pairs. 

The pairs gyrating in a transversal magnetic field have a flat 
distribution in 8 (not cos 8). Therefore, their synchrotron emission 
has a pattern 

^— = * i-U (23) 
d cos 8' 4-7T 7r sin 8' 

An identical distribution would be produced by SSC emission (be- 
cause the synchrotron photons that are Comptonized by 2D pairs 
are mostly produced by the next generation, lower energy, isotropic 
pairs). With the relation sin 8' = 8 sin 8, we get the amplification 
factor (see dashed curves in Fig. [8jl 



(24) 

This distribution has sharp peaks at 8 — and n. 

The external radiation IC scattered on 2D pairs produces a 
pattern 

AL' _ L' 2 (l + cos(9') 2 
d cos 8' 4-7T 3ir sin 8' ' 

which gives the distribution of the observed emission with the 
peaks at 8 — 0, n (see dotted curve in Fig.[8]l: 

2 S 4 1 



MP) 



3tt Tf sin 8 ' 



(26) 



3.6.3 Radiation from external environment 

In addition to the radiation from the jet, there is also radiation from 
the relativistic pairs produced in the external environment. They 
also have two populations. The first population is produced by the 
high-energy photons from the jet propagating a small angles to the 
jet direction 9 < 1/T. They perform 2D circular motion in the 
transversal magnetic field of the surrounding. The second popula- 
tion produced by the pair cascade is nearly isotropic. 

The synchrotron radiation of 2D pairs follows distribution 
J23t , with 8' replaced by 8 (see short-long dashes in Fig.[8). 

(27) 

7r sm8 

The IC scattered isotropic BLR photons have the same distribution. 

However, because most of the seed photons (synchrotron pho- 
tons from the jet or the accretion disc photons) are beamed in the 
direction of the jet, the IC emission by the external pairs follows 
the distribution l |25t with 8' replaced by tv — 8 (see long dashes in 
Fig.[8}: 

2 (1 -cosfl) 2 



L{8) oc 



3tt 



(28) 



Thus, the 2D pairs produce a sharp peak at 8 — tt. 

The emission of isotropic pairs Comptonizing the disc (or jet 
synchrotron) radiation follows the law J2U with 8' replaced by n — 
8 (see dot-dashed curve in Fig. [8}: 



L{8) cc |(1 -cos(9) 2 . 



(29) 



And, finally, isotropic pairs scattering isotropic radiation obviously 
produce a flat distribution of L(8). 

3.6.4 Emission pattern: simulations 

As an example of the angular distribution of radiation predicted by 
the photon breeding model, we present results for run 01. The his- 
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Figure 9. Electron (and positron) energy distributions in the jet comoving 
frame averaged over the jet volume for runs 01-04 and 10-15. 



tograms in Fig. [8] show the dependence of the luminosity in three 
typical energy bands (IR-optical, X-ray and gamma- rays). We see 
that the distribution of high-energy photons at small 8 follows quite 
closely that expected from the jet of Tj = 20 emitting isotropically 
in its rest frame (e.g. SSC, given by equation |20ll. The ERC radia- 
tion from the jet is more beamed (see equation [22] and the red solid 
curve in Fig.[8j. An excess at 8 — can be explained by the contri- 
bution of the SSC and ERC emission from 2D pairs (given by eqs. 
|24]and|26] respectively; see dashed and dotted curve in Fig.[8j. 

At angles 8 > 0.45, another component is obvious. This large- 
angle component is nearly isotropic and its flux at infinity is 4.5 or- 
ders of magnitude smaller than the average flux observed within the 
jet opening. For the run 01 with Tj = 20, the observed luminosity 
of this large-angle component at 8 = 7r/2 exceeds the luminosity 
given by expression j20l by more than 3 orders of magnitude. This 
off-axis emission is produced by the pairs gyrating in the external 
environment via Compton scattering of the disc and BLR photons 
and synchrotron photons from the jet. A peak at 8 — tt can be well 
described by IC emission of 2D pairs in the external environment 
(see long dashes and short-long dashes in Fig. [8} given by formulae 
<27> and J28>. 

The low-energy (IR-optical and the X-rays) photons are much 
less beamed. The spatial gradients of V is the main cause of that. 
For example, the average V of the inner half of the jet volume with 
r < 0.7 in run 01 is about 18, while the outer 20 per cent move 
with F < 10, dropping to F ~ 2 in the outer 5 per cent (see also 
Figs|4]and[5]l. 

The soft band is dominated by the synchrotron radiation 
(Fig.[7). The emission at 8 = 0.5 exceeds the simple estimate d201 > 
by 3 orders of magnitude. At angles 8 > 1/Fj = 0.05 the emission 
is mostly produced in the slower, decelerated layers of the jet (see 
Fig. |4), and the emission at 8 > 0.5 is produced in the layers of 
r ~ 2. The synchrotron from 2D pairs in the jet even produces a 
peak at 8 — n. Some contribution to this peak is provided by the 
IC (low-energy tail) emission from the external medium. 

The emission in the X-rays is the mixture of the synchrotron 
and ERC emission resulting in the intermediate behaviour seen in 
Fig. [8] The large-angle emission here is produced by IC in the ex- 
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Figure 10. Kinetic luminosity carried by pairs with respect to the initial jet 
power for runs 21-24 (Tj =40-14). 

ternal environment (see Fig.0 and therefore has an angular depen- 
dence similar to that of the gamma-rays. 

The sharp feature at 9 = will be washed out if the jet is 
actually conical rather than cylindrical or the jet magnetic field is 
tangled. The tangled magnetic field in the external medium will 
also smooth out the predicted peak at 6 — n. 

3.7 Pairs produced in the jet 

The electron and positron comoving energy distributions averaged 
over time and jet volume for runs 01-04 and 10-15 are shown in 
Fig. [9] All spectra except that for run 15 have a clear low-energy 
cutoff at 7 C = 10-10 3 , this is the energy the pairs cool down to 
in the dynamical time. The cutoff comoving energy assuming the 
Thomson regime for Compton cooling and neglecting SSC can be 
expressed as 



,c t dyn (£ B +£ lso r 2 y 

where idyn ~ 10/r ~ 1 is a typical (comoving) time the particles 
spend in the active cylinder. In our examples the Compton losses 
dominate over the synchrotron losses (see equation [9] and Table 1). 
For runs 04, 14 and 15 the cutoff is smoother because the cooling 
rate above 7 C ~ 3 x 10 3 is affected by the Klein-Nishina regime, 
making it slower. 

The spectrum above the cutoff is the standard cooling spec- 
trum diV/d7 oc 7~ 2 , and it changes to a steeper slope at higher 
energies. This is the result of the pair cascade developing in the jet. 
The cascade pairs affec t the energy dis tribution making it steeper 
dN/d'y oc 7~ 3 (e.g. ISvenssoal 19871) above the break energy 
7b ~ 10 -10 . Pairs of lower energies produce photons that can 
freely escape from the jet and the cascade stops. We may note here 
that most of the pairs in the jet are produced by photons emitted 
in the jet. Such process does not contribute to the energy release. 
However, during the cascade the typical photon energy decreases; 
it leads to the increase of the probability of the photon escape to the 
external environment. This in turn increases the efficiency of the jet 
energy dissipation. 
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Figure 11. The spectra of the external radiation and of the soft (mostly syn- 
chrotron) radiation induced by the jet activity for run 3 1 (solid curves) and 
run 35 (dashed curves). Notice that the spectra here are E AN/AE, not 
E' 2 AN/AE, to emphasise the number density of target soft photons avail- 
able for pair production. Upper panel: the collimated components, lower 
panel the isotropic external radiation and the large-angle (8 > 40°) syn- 
chrotron jet component. 

Fig. [To] shows the fraction of the jet energy carried by 
generated pairs, which reaches 2-5 per cent. The number ratio 
of generated electrons and positrons to that of the protons in- 
creases with the jet power and for runs 14-10 it is t±/t p — 
0.017, 0.135, 0.63, 1.9, 4.0, respectively. 



4 PROPERTIES AND LIMITATIONS OF THE MODEL 
4.1 Non-linear effects 

There are two kinds of non-linearity in the problem: one appears 
because of the feedback of radiation on the fluid, another is the re- 
sult of the synchrotron emission, which produces additional opacity 
for the high-energy photons. Here we consider the second one. The 
energy density of soft synchrotron radiation of pairs produced in 
the outer decelerated layers of the jet can exceed that of the exter- 
nal radiation described in Section l2.3.2l In Fig.QTTwe compare the 
primary external radiation to the jet emission averaged over the vol- 
ume. For run 31 the latter significantly exceeds the external radia- 
tion in both components: longitudinal and isotropic. In this case the 
synchrotron emission of the jet becomes the main target for high- 
energy photons. For runs 01, 02, 10-12 and 35, the soft radiation 
associated with the jet is comparable with the original external radi- 
ation. In runs 03, 13-15 (also runs 41 and 42 for model B), despite 
rjj being the same, the jet emission at low energies is much smaller 
than the external radiation because the jet spectrum is harder. These 
cases can be treated as linear, so that the solution can be scaled to 
any lower jet power. 

In the case of substantial non-linearity, the process becomes 
self-supporting. The external radiation becomes unnecessary, be- 
cause the stronger sy nchrotron radiation plays the same role. 
IStawarz & Kirkl d2007t) have studied this kind of non-linearity us- 
ing kinetic equations and found that a source of high-energy pho- 



High-energy emission of relativistic jets 13 



tons quenches itself at high compactness. In our case, we see that 
the efficiency decreases at high compactness and the final spectrum 
gets softer. We can expect that such a tendency will continue with 
the increasing compactness. 

The synchrotron emission from the jet scattered in the BLR 
(neglected i n our simulations) can als o affect the development of 
the cascade dGhise llini & Madau 1996). The dominant contribution 
to the scattering is provided by the region of thickness ~ Rj imme- 
diately adjacent to the jet. If we associate the fraction of radiation 
scattered/reprocessed in BLR, 771 , with the optical depth through 
the region of size R, then the optical depth through the layer of 
thickness Rj for synchrotron photons moving at angle 1/Yj is of 
the order T sc ^rii9Yj. We can write now the scattered luminosity as 



to the kinetic energy of the flow one gets the final velocity of the 
entrained medium (in non-relativistic limit): 



i s , blr = L b t bc = ■q s {-q ] L d )r) i 9T i , 



(3D 



where r/ s = L s /Lj < 0.1 is the fraction of the jet power emitted 
in the form of synchrotron radiation. Because the scattered syn- 
chrotron is escaping through the surface of ~ A-k8 R 2 , the ratio of 
its energy density to that of the BLR photons is 



f/s,BLR 

Ublk 



(32) 



This shows that the scattered synchrotron can be important for large 
m orI Y 



4.2 Temporal variability 

In our simulations, we have assumed that the jet power and the 
Lorentz factor at the inlet do not change. However, because of the 
large non-linearity, the output photon flux demonstrates large non- 
Poisson fluctuations at different time-scales. At the present level 
of numerical resolution, it is difficult to distinguish confidently be- 
tween numerical fluctuations and natural chaotic behaviour of the 
non-linear dynamical system. We just would like to emphasise that 
the latter possibility is very probable. 

The transversal gradient of the Lorentz factor, dY/dr, affects 
the rate of the photon breeding in the same way as the slope of a 
sand pile affects sand avalanches. The emitted radiation feeds back 
on the jet dynamics, making the gradient smoother. Thus dY/dr 
balances near its critical value in the same way as the slope of the 
sand pile does, when one pours sand on the top. It seems that we 
observe this in the cases of highest non-linearity, when the soft ra- 
diation of the jet dominates the external soft radiation (runs 10 and 
31, see Fig. II It. It these cases the radiation power is more sensi- 
tive to the transversal gradient of F because of a short free path of 
high-energy photons in this direction. 

If we assume a large fluctuation in the jet power at a time- 
scale shorter than R/cTf, the photon avalanche will be propagating 
along the jet and the observer will see a large increase in the lumi- 
nosity partially because of the time compression. Simulations of 
strongly variable sources in the context of photon breeding model 
are left for future studies (see also SP06). 



4.3 Impact on the external environment 

The high-energy photons produced in the jet interacting outside 
of the jet deposit their momentum into the external environment 
and therefore accelerate it. This acceleration can be neglected if 
the density of the external matter is sufficiently high. Let the mo- 
mentum deposited in a unit volume per unit time is k(t, z) = 
dp/(dVdi). Then, equating the work of the photon pressure force 



P(r) = v/c 



1/2 



(33) 



where n cm is the particle density of the external medium. The 
required condition /3 <C 1, allows to put the lower limit of the 
density of the environment: 



t! » tier = 



2 J K(r,z)dR 



(34) 



During the simulations we accumulate the average deposited 
momentum k. Its dependence on r is strong. For the layer closest to 
the jet, 1 < r < 1.005, we have / «(1 < r < 1.005, z)dR = 17 
erg cm for run 31, which translates to 10 cm . The 

requirement seems reasonab le for the BLR, where the cloud 
density is about 10 cm" ~ 3 dOsterbrock & Ferland l2006) and even 
the mean density should not be much smaller. The critical densities 
for other runs are even smaller: n cr = 10 3 cm -3 for runs 02 and 
12, n C r = 30 cm -3 for runs 04 and 14. For runs 41 and 42, the 
required density is n cr = 10 cm -3 , which is still reasonable at a 
parsec scale. 



4.4 Important phenomena beyond ballistic approximation 

Our consideration is exact and model independent when concerning 
photons and pairs: their interactions, propagation, breeding, etc. at a 
given distribution of the fluid Lorentz factor and a given magnetic 
field. The only assumption we make is the spectrum and density 
of seed high-energy photons, this is however not important as their 
energy content is negligible and the system "forgets" their spectrum 
very soon. 

On the other hand, our consideration is very approximate and 
is based on a simplified model when it concerns the fluid and the 
feedback of the particles on the fluid. How our simplifications af- 
fect the results and what we would obtain with a more realistic 
hydrodynamical treatment of the problem? 

First, the pressure will redistribute the Lorentz factor of the 
fluid. This will hardly change significantly the process of particle 
breeding. More important, the formation of internal shocks at re- 
gions of maximal deceleration can increase the radiation efficiency 
by accelerating charged particles. 

Then, we can expect that a turbulent layer near the jet bound- 
ary will be formed due to the large shear gradient of the Lorentz 
factor. The same can happen due to the Kelvin-He lmholtz instabil- 
ity of the boundary discontinuity l lDe Y oung 1986), however for an 
ultra-relativistic jet this instability develops slowly. The radiative 
deceleration of the outer layer of the jet can generate the turbulent 
mixing layer much earlier. Then the turbulent layer should again 
contribute to the radiation d ue to diffusive particle acceleration (see 
IStawarz & Ostrowskill2003h . 

Finally, we ignore adiabatic heating of relativistic electrons 
while the jet (comoving) compression ratio in z direction is ~ 
Tj/r, where Y is the Lorentz factor of the decelerated fluid. The 
account of the adiabatic heating is not simple, because the jet can 
expand in the transversal direction. The compression also leads to 
the amplification of the magnetic field Bj oc 1/T (ignored in our 
simulations), which would in turn increase the role of synchrotron 
radiation in the decelerated part of the jet. 

The formation of shocks, turbulence and adiabatic heating 
of electrons, which can be caused by the jet deceleration, should 
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Figure 12. Fraction of the jet power converted into photons versus time for 
run 30. The initial density of the high-energy, seed photons corresponds to 
the extragalactic gamma-ray background. 

increase the radiative efficiency of photon breeding mechanism, 
while the magnetic field amplification probably decreases it. These 
phenomena are interesting and should be studied in future simula- 
tions. 



4.5 Origin of the seed high-energy photons 

In the presented simulations the energy content of the seed high- 
energy photons, which initiate the further evolution of the system, 
was taken to be 5-7 orders of magnitude smaller that the final 
steady-state total energy content in photons within the jet. Where 
these initial photons can originate from? They can be associated 
with the shock acceleration of charged particles in the jet. There 
is, however, a minimum possible source of high-energy, seed pho- 
tons provided by the extragalactic gamma-ray background. The ob- 
served gamma-ray backgroun d intensity at Earth at e > 10 is 
~ 10" 6 GeV cm" 2 s" 1 sr" 1 dSreekumar et al.lll998l) . correspond- 
ing to the energy density 6 x 10~ 19 erg cm -3 . This would give 
the energy ~3x 10 31 i?? 16 erg within the considered cylindrical 
volume of radius Rj and height 20Rj . 

We thus made a simulation (run 30) where the total energy 
of the seed high-energy photons in the simulation volume was 20 
orders of magnitude smaller (approximately corresponding to the 
energy of the extragalactic gamma-ray background) than the jet en- 
ergy in the same volume Ej ~ 10 51 erg. The result is shown in 
Fig. [12] The system passes 20 orders of magnitude in the total pho- 
ton energy during 300iij /c ~3 years. A jump at the very begin- 
ning of simulations corresponds to almost immediate conversion of 
seed photons into pairs with the corresponding energy increase by 
a factor rf = 400. The photon energy grows exponentially with 
the time constant t p b = 2.7 at t < 25 and t p b = 8.2 at later 
time. The rapid rise at t < 20 corresponds to the growth of the 
photon avalanche as it moves downstream with the flow. At later 
times the growth is supported by the spatial feedback loop due to 
up-streaming photons. The final steady-state regime does not de- 
pend on the initial density of high-energy photons. This numerical 
experiment means that a jet, with parameters corresponding to run 



30, would radiate half of its total energy in three years just due to 
the photon breeding starting off the extragalactic gamma-ray back- 
ground. 



5 ASTROPHYSICAL IMPLICATIONS 

5.1 Spectral energy distribution of blazars 

Blazars demonstrate two-component spectra. The low-energy com- 
ponent peaking in the IR-X-ray range (the flares in the TeV blazars 
show hard X-ray spectra peaking above 100 keV) is tradition- 
ally interpreted as synchrotron radiation from relativistic electrons. 
The high-energy component is the Comptonized radiation from 
some seed soft photons. The spectra of low-luminosity objects 
(TeV blazars, BL Lacs) are well described by the synchrotron 
self-Compton (SSC) models. In high-luminosity blazars, both syn- 
chrotron and Compton components peak at lowe r energies (so 
called blazar sequence, iFossati et al. I I 19981 ; but see iGiomrm et al.l 
l2007h . The role of the external radiation (from BLR or dusty torus) 
in providing seed photons is probably more important, and the 
spectra can be described by ex ternal radiation Compton (ERC) 
model dSikora et al.lll994 fl997). In all objects, the two compo- 
nents are rather well separated. 

The spectra, predicted by the photon breeding model (see 
Fig. [6^), demonstrate similar features. At low luminosities and 
compactnesses (runs 04, 14-15), the mean Lorentz factor of the 
electron-positron distribution is large (see Fig. |9). In spite of the 
fact that (in the jet frame) the magnetic energy density is 1-2 or- 
ders of magnitude smaller than that of the isotropic radiation (i.e. 
(■b <C fi sc ,r 2 , see Table 1), the synchrotron power strongly domi- 
nates over the external Compton, because the external photons in- 
teract with the high-energy pairs in the Klein-Nishina regime. As 
most of the synchrotron photons are also above the Klein-Nishina 
cutoff ~ 1 /y, the SSC is weaker than the synchrotron. 

At higher luminosities (and compactnesses), the high-energy 
photons produce a pair cascade. This reduces the mean energy per 
pair (see Fig. [9] runs 01-03, 10-13). At these compactnesses the 
pairs also cool more efficiently producing a cooling distribution 
cL/V/dy oc 7~ 2 (see Section [3~7t , which changes to a steeper 7~ 3 
behaviour at higher energies as a result of pair cascade. The latter 
part of the pair energy distribution is responsible for the flat part of 
photon spectrum (in E 2 AN/AE units)[jThe reduction of the mean 
pair energy results in a shift of the photon spectra to lower energies 
(Fig.|6ji). It also leads to the increased role of the ERC, because now 
most of external photons interact with pairs in Thomson regime. 

The properties of our model are the increasing role of the 
ERC over SSC and synchrotron with luminosity and a shift of the 
spectral peaks to lower energies. These are in general agreement 
with the observational trends. However, there are some discrepan- 
cies between the shapes of the theoretical and the observed spectral 
shapes. At low luminosities, the synchrotron peak in our model is in 
the 100 MeV-1 GeV range, while it is normally observed in the X- 
ray range. One need to note here that the angle-dependent spectra 
at 6 > 1/Fj have sharper features (see Fig.[7J and the synchrotron 
spectrum peaks at slightly lower energies that the spectrum aver- 
aged over 10° around the jet axis. Another possible reason for the 
discrepancy is that the electron distribution is broad and extends 

5 One should note here that the pair cascade does not degrade the pair 
energies infinitely, when the energy of the photon in the jet frame drops 
below l/SmaxTj ~ 10 4 , it can easily escape and the cascade stops. 
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to very high energies. This is a consequence of our assumed spec- 
trum of the soft isotropic radiation extending to a very low energy 
x m in = 1CF 9 . A narrower and harder distribution of external soft 
photons might make the injection function of pairs narrower and 
softer, so that the maximum pair energy is below 7 max given by 
equation dl3t , reducing thus the energy of the synchrotron cutoff. 

At high luminosities, our model predict basically a single 
component where the synchrotron and ERC components are mixed, 
while there are clearly two separated components in the data. How- 
ever, there is no consensus whether the two components are co- 
spatial. For example, the bright blazar 3C279 did not show al- 
most any variability in the radio-optical band, when the gam ma-ray 
flux changed by an order of mag nitude dWehrle et al] [l998). Thus 
the radio-to-optical emission in the high-luminosity sources can be 
produced is a separate region, far away from the black hole. 

In our model, the pairs escaping from the region filled with 
soft photons still carry a considerable (a few per cent) fraction of 
the jet energy (see Fig.UOt. If these pair are reheated at larger scales 
(e.g. by a diffusive acceleration by various plasma perturbations), 
where soft radiation field is weaker, they will radiate most of their 
energy in the form of synchrotron. At the expense of 3-5 per cent 
of the jet energy for pair reheating, the synchrotron peak observed 
in low-peak blazars could be reproduced. If this hypothesis is valid 
then the high- and the low-energy peaks have different origin: the 
high-energy peak is a mixture of prompt synchrotron and Comp- 
ton radiation of high-energy pairs in the fast cooling regime. The 
low-energy peak is the synchrotron radiation of pairs in a heat- 
ing/cooling balance. In such a scenario the intensity of high- and 
low-energy peaks should correlate in time, but in a loose manner, 
the correlation should be stronger at longer time-scales. 

5.2 Electron energy distribution in BL Lac objects and 
particle acceleration 

One of the most often used in the recent years scenario for for- 
mation of the relativistic electrons, responsible for the observed 
emission, involves the diffusive (Fermi) acceleration at relativistic 
collisionless shocks (or shear layers). It has been claimed that the 
acceleration at relativistic shocks produces the power-law particle 
spectrum dN/dq oc q~ p with a universal p = 2 .2 index (see e.g. 
lAchterberg et al.ll200ll ; lKeshet & Waxmanll2005l) which is close to 
the v alue observed in BL Lacs (Ghisellini . Celotti & Costamantd 
2002). These results, however, were based on certain simplify- 
ing assumption regarding the scattering of particles by magnetic 
inhomogeneities. Using accurate Monte C arlo simulations with 
a mor e realistic magnetic field structure, Niemiec & Ostrowskil 
d2006h demonstrated that no universal power-law is produced and 
questioned the role of the Fermi-type particle acceleration mecha- 
nism in relativistic shocks. 

Using particle-in-cell simulations, ISpitkovskvl i f2007t) recently 
showed that the electrons can be heated by interaction with ion 
current filaments in the upstream region reaching almost an en- 
ergy equipartition with protons. The low-energy cutoff appears at 
7min ~ irm p /m c . In internal shock model, the relative veloc- 
ity of colliding shells is mildly relativistic and thus we expect 
7min ~ 500. At this stage, it is too early to say whether the model 
would be able to produce a power-law electron distribution extend- 
ing to very high energies. 

The detailed fitting of the broad-band spectra of TeV blazars 
and low-power BL Lacs with the SSC model shows that the distri- 
bution of relativistic electrons injected to the system can be mod- 
elled as a power-law of index p « 2-2.5 with the low-energy cutoff 



at 7 m in ~ 10 4 -10 5 an d the high-energy cutoff at 7 max ~ 10 6 -10 7 
dGhisellini et alj|2002t iKrawczvnski et al J [20021 ; iKonopelko et al.l 
l2003l;lGiebels et alj|2007l) l 6 l Such a peaked distribution of injected 
electrons is not consistent with a diffusive process. 

The photon breeding mechanism, on the other hand, predicts 
the injection spectrum to be bounded and to mirror (relative to 
m c c 2 ) the spectrum of the soft photon background, because the 
pairs in the jet are produced only by those high-energy photons 
that can interact with the soft photons. For example, if soft pho- 
tons are from the accretion disc of maximum temperature max , 
the low-energy injection cutoff is at 7 m i n ~ Vj/6Q ma x (see SP06; 
factor Tj comes from the Lorentz transformation to the jet frame). 
This cutoff should depend weakly on the luminosity and the black 
hole mass, which defines the characteristic emission radius, so that 
Omax » 10" 5 Ly4 5 (MBH/10 8 M Q )" 1/2 and thus we can predict 
that 

At high compactnesses, however, the synchrotron photons from the 
jet can provide enough opacity (see Section |4~T| | and 7 m in may be 
smaller. 

For the conditions appropriate for AGNs, the high-energy 
cutoff of the pair distribution 7 ma x, that can participate in pho- 
ton breeding, is defined by the different physics. An electron (or 
positron) of Lorentz factor higher than that given by equation d 1 3 1 > 
loses most of its energy to the synchrotron emission which escapes 
freely from the jet. Thus, at low compactnesses the pairs are in- 
jected to the jet between 7 m i n and 7 ma x. At high compactnesses, 
because of the higher opacity the Compton scatted photons cannot 
easily escape from the jet, producing a pair cascade which degrade 
their energy (see Fig. |9). The range of Lorentz factors of injected 
pairs as predicted by the photon breeding mechanism is in general 
agreement with observations of BL Lacs. 

5.3 Jet power and composition 

In many models of the high-energy emission from relativistic jets, 
the radiative efficiency is low because only internal energy of the 
jet can be dissipated (e.g. in internal shock model) and only a small 
fraction of the energy is transferred to the (high-energy) electrons. 
This leads to the large jet power requirements. 

The photon breeding mechanism has important properties that 
significantly reduce the energy budget required: (i) it produces only 
electrons (positrons) of high Lorentz factors; (ii) it efficiently de- 
celerates the jet, using thus most of its total power. Thus the radia- 
tive efficiency is extremely high, reaching in many cases 50-80 per 
cent. 

The primary composition of the jet is not constrained if the 
photon breeding works, except the general requirement that it 
should not be dominated by cool pairs. Indeed, a dynamically im- 
portant component of cold pairs can be excluded because of the 
absence of the bulk C ompton scattering bump p roduced by such 
pairs in the soft X-rays dSikora & M adeiski 2000). Moreover, large 
contribution of pairs to the bulk energy might be problematic, be- 
cause of their short annihilation time in the vicinity of the black 

6 The high-luminosit y objects require much lower energy electrons 
iGhisellini et"ai]|2002h . but this result may be biased because the electron 
distribution is obtained from one-zone models, while the synchrotron peak 
in reality may be not related to the gamma-rays. 
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hole {piandford & Levinsorill 19951) . Magnetically dominated jet is 
less favourable for photon breeding at V < 20, while for a larger 
Lorentz factor the radiative efficiency is not sensitive to the jet mag- 
netisation. 

Independently of the primary jet composition the photon 
breeding loads the jet with a significant amount of relativistic pairs. 
The contribution of these pairs to the jet energy flux according to 
our simulations is 2-5 per cent (see Section 13771 . This is probably 
an underestimation as we do not take into account diffusive acceler- 
ation and adiabatic heating of pairs. The number ratio of generated 
electrons and positrons to protons (in the case of matter-dominated 
jet) may be related to the mean electron random Lorentz factor (as 
measured in the jet frame) and the radiative efficiency: 
r± m p £ off 

~ TT- (-to) 

r p m c (7) 

For low-power jets, a typical 7 ~ 10 5 -10 6 , and the number of gen- 
erated pairs should be relatively small. At high powers, the mean 
pair energy decreases and the generated pairs dominate (Section 
13. 1\ . Thus, for high-power quasars the model predicts the number 
of pairs exceeding in number the protons, while dynamically the 
protons still dominate. This is consistent with the constraints on the 
pair content obtained from the X-ray s pectra of blazars associated 
with quasars dSikora & Madejski 200(J. 



ing stops. If at the parsec scale, the jet is still relativistic enough, the 
process can start operating again on the IR radiation from the dusty 
torus. The typical dust temperature is Gdust = fcXkust /m c c 2 w 



-1/4 p-l/2 
J d,45 JX P C 



The pair-production optical depth is about 

Ld,45 



= n ph , dus t(0.2<T T )flj = 30^(200) 



-^pc@dust,- 



30^(208) L 3 d %R; 



-1/2 



(37) 



where n p h,dust =0.37 t/dust / fcTdust is the number density of pho- 
tons from the dust and rji is the ratio of the dust to the disc lumi- 
nosities. Thus, even for r\\ of the order of a few percent, the photon 
breeding is effective at the parsec scale in bright quasars. 

It has not escaped our attention that the jet may become active 
even further out: by interacting at kpc scale with the stellar radia- 
tion field and at ~ 100 kpc with the cosmic microwave background 
radiation. 

We note here that photon breeding works well at high disc lu- 
minosities (see Section [33) characteristic for quasars. In BL Lac 
objects, where the discs are less luminous, photon breeding may 
operate at small distance scale, which in turn requires small black 
hole mass and/or fast jet acceleration, or for large jet Lorentz fac- 
tor r, > 40. The obs erved fast vari ability from PKS 2155-304 
dAharonian et alj|2007h and Mrk 501 jAlbert et al]|2007h supports 
these suggestions. 



5.4 Conditions for photon breeding and gamma-ray emission 
sites 

The arguments based on the gamma-ray transparency and the ob- 
served short time-scale variability constrain the gam ma-ray emit- 
ting r egion in blazars to distances about ~ 10 17 cm dSikora et al.l 
1 1994b . This distance is often associated with the BLR. In the in- 
ternal shocks models, the shells of Lorentz factors Ti and Y2 
ejected from the central source and initially separated by the dis- 
tance R collide at the distance R = RoVi/[l - (Ti/r^) 2 ]. Tak- 
ing R = 10 15 cm (« 30i? s for a 1O 8 M black hole), I\ = 10 
and T2 3> Ti we get R — 10 17 cm. However, the agreement with 
the aforementioned constraints is rather coincidental, because the 
black hole mass can vary by 3-4 orders of magnitude and there is 
no physical theory yet predicting the assumed separation between 
the shells. 

The photon breeding mechanism allows various emission 
sites. If the jet is already accelerated and collimated at distances 
~ 100/?s from the black hole, the direct radiation from the ac- 
cretion disc can serve as a target for high-energy photons for pair 
production and to trigger the photon avalanche. The photon-photon 
opacity through the jet is very large at these distances and the jet 
deceleration and radiative efficiency might be rather small, because 
only a very thin boundary layer is active. The escaping radiation 
should have a high-energy cutoff at tens of GeV rather than at 
TeV energies. At larger distances the disc radiation becomes more 
beamed along the jet and the disc photons do not interact any- 
more with the high-energy radiation produced in the jet. Thus the 
photon breeding stops. The jet now is, however, loaded with high- 
energy photons that can serve as seed high-energy radiation for the 
avalanche to be triggered at larger distances, where the conditions 
for photon breeding are satisfied. The natural scale for that is the 
position of the BLR, where the density of the isotropic soft pho- 
tons is sufficiently high. As we showed in Section [3] the radiative 
efficiency here can be very large. 

At even larger distances, the BLR photon density drops, the jet 
becomes transparent for high-energy photons and the photon bread- 



5.5 The Doppler factor crisis in TeV blazars 

The Doppler factors required by the homogeneous SSC models 
to describe the spectra of the blazars emitting at TeV energies 
(to avoid absorption b y the IR radiation), is rath e r large (5^20-50 
1 Ghisellini et alj|2002l: Krawczvnski et al J |20o3 : iKonopelko et al.l 
l2003l : lGiebelsetal.ll2007h . Such high 5 values disagree with the 
small apparent velocities observed at the parsec scale i n Mrk 421 
and Mrk 50 1 jMarscheJll999[: IPiner & Edwards] l2005h and other 

TeV blazars (Piner & Edwards 2004). 

One proposed solution to this 5-crisis dTavecchiol 120051) is 
that the jet Lorentz factor drops from 15 at the gamma-ray emit- 
ting, subparsec scale to ~4 at the ra dio-emitting VLBI parsec scale 
( Georganopoulos & Kazanas 2003). For the viewing angle within 
the beaming pattern of the initial jet, 9 < 1/Ij, the apparent veloc- 
ity of the decelerated jet with Lorentz factor Tf <C Tj is 



Pi sin 6 
1 — Pi cos I 



26>rf 



1 + (evtY 



r 2 

26T? < 2— 



(38) 



Thus, deceleration of the jet to Tf < ^Y- s /2 guaranties that the 
apparent motion is subluminal. 

The question arises what is the mechanism that de celerates the 
jet at subparsec scale. The internal shocks in the jet 1 Spa da et al.l 
l200lh are rather inef ficient unless huge fluctuations of the Lorentz 
factors are involved teeloborodovl2000l) . which seems unlikely. On 
the other hand, the photon breeding mechanism predict significant 
deceleration of the jet (see Table 1) and strong emission at subpar- 
sec scale, if the conditions for photon breeding are satisfied. The 
photon breeding thus may be responsible for jet deceleration and 
resolution of the 5-crisis. 

A different way to solve the 5-crisis i s to assume that the jet 
consists of a fast spi ne and a slow sheath (Chiaberge et alj|2000l : 
Ghisell ini et al . 2005). In this model a slower sheath having a softer 
spectrum dominates the emission in the radio band. The photon 
breeding mechanism produces exactly such a structure (Fig.|4j and 
the spectra (Fig. [8} in a self-consistent way. 
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Thus, the jet behaviour predicted by the photon breeding 
mechanism unites the two proposals: it decelerates the jet at sub- 
parsec scales and produces a fast spine-slow sheath structure. 



5.6 Emission pattern 

5.6.1 Unification of BL Lacs and FR I radio galaxies 

Acc ording to the unificatio n scheme of radio-loud active galax- 
ies dUrrv&Padovanill 19951) . BL Lac objects are Fanaroff-Riley I 
(FR I) radio galaxies with their jets oriented close to the line of 
sight. The average Lorentz factor of the jets, derived by match- 
ing the luminosity functions and statistics of BL Lac and FR I 
(and B 2) samples, was estima ted to be Tj ~3-5 dUrrv & Padovanil 
1 19951 : lHardcastleetai]|2003l) . The core optical luminosities of 
BL Lacs are on avera ge 10 4 ' 5 times higher than those of FR I 
dChiaberge etai1l200Cl) with the same extended radio luminosity. 
Associating the core emission with the jet, we again get rather low 
Lorentz factors Tj = 6 or 4 for the continuous steady jet or a mov- 
ing blob, respectively. These low Fj are in clear disagreement with 
the much higher values of the Doppler factors required to explain 

the gamma-ray emission of BL Lacs. 

Both decelerating ( Georganopoulos & Kazanas ]2003h and 
structured dChiaberge et alJl200Ct iGhisellini et alj|2005h jet mod- 
els provide explanations for the discrepancy. Because the photon 
breading mechanism produces decelerating and structured jet, it ob- 
viously is also capable of unifying BL Lacs with the radio galaxies, 
keeping high Tj to explain the gamma-ray emission. In blazars the 
emission from the fast, Fj ~ 20, spine dominates. In radio galaxies 
observed at large angles 9 ~ 1 (see Fig.(8](, the optical emission is 
produced in the slower sheath of F —2-6. 



5.6.2 Off -axis emission and the TeV emission from radio galaxies 

As we have discussed above, the photon breeding mecha- 
nism naturally produces very broad beams of photons (see als o 
iDerishev et alj 120031 ; iDerishev. Aharon ian & Ko charovskvl |2007|) . 
The gamma-ray emission at 6 ~ 1 is produced in the external en- 
vironment (see Fig. [8). The luminosity ratio between the nearly 
isotropic emission and the beamed emission at 6 w 1/Tj is about 
ry 4 (two powers of Fj come from the energy amplification in 
the jet and two powers appear because of beaming). This ratio is 
~ r?/10 larger than that predicted by equation l |20| >. Thanks to 
the emission from the external medium, the nearby misaligned jets 
become observable in high-energy gamma-rays. 

The central source of M87, which has the best studied 
jet, has a rather low d i sc lum inosity of only ~ 10 42 erg s _1 
teiretta. Stern & Harris! Il99ll) . With the threshold L d > 
10 43 i?i7erg s" 1 for T = 20 the photon breeding in model A is 
possible at R < 10 16 cm, which is only ~ 10i?s for the estimated 
mass of the central black hole ~ 3 x 1O 9 M0 dMacchetto et al.1 
1 19971) . At such a distance the photon breeding will operate on the 
direct disc radiation (model C). The implied distance is consistent 
with the detection of th e rapid (~ 10 s) var iability of the TeV 
photon flux from M87 (Aharonian et al. 2006). However, photon 
breeding may not ope rate, if the wide p attern a t the base of the 
jet ob served by VLBI dJunor et alii 1 999l but see lKrichbaum et ail 
1 20071) actually reflects the jet geometry, which would mean that the 
opening angle is large. 



5.7 Observational appearance and volume dissipation 

Most of the emission from the jet undergoing photon breeding 
comes from the regions of maximum gradient of F. In the steady- 
state, the volume emissivity contributing to the observed emission 
(at small angles to the jet) is 



j(r, z) oc 



dF(r, z) S 3 
dz ~T 



(39) 



In the case of small luminosities (and small Fj) only the outer lay- 
ers of the jet suffer significant deceleration (see Fig.|4jj, run 24), re- 
sulting in the limb-brightening, with most of the emission coming 
from large z, where the cascade had time to develop. The situation 
is more complicated at high luminosities (or large Tj). The outer 
layers decelerate rapidly (see curves for runs 21-23 in Fig.|5j mak- 
ing the emission limb-brightened at small z, while the layers close 
to the jet axis decelerate slower. Thus the emission from the core 
dominates at large z. At angles larger than > 1/Tj, only slower, 
significantly decelerated layers contribute to the observed luminos- 
ity, and therefore we always would see the limb-brightened emis- 
sion. 

Our predictions are consistent with the observed limb- 
brightened morpholo gy of the Mrk 501 jet at the parsec scale 
dGiroletti et"ai]|2004l) . Such a structure could a result from the ef- 
ficient deceleration and loading by relativistic pairs of the jet outer 
layers by the photon breeding at a subparsec scale. 



6 CONCLUSIONS 

The photon breeding mechanism can be an extremely efficient way 
of dissipating bulk energy of relativistic (Fj > 10) jet into high- 
energy radiation. In the case of high luminosity AGNs and power- 
ful jets, the mechanism is very robust: it works at any value of the 
jet magnetic field and is not sensitive to the spectrum of the external 
radiation. At high luminosities and jet Lorentz factors, the mecha- 
nism should work almost inevitably, the only fact that can prevent 
the runaway photon breeding is a purely longitudinal (along the jet) 
geometry of the magnetic field. 

At the intermediate range of AGN luminosity, 10 43 < 
LdRn < 3 x 10 44 erg s _1 and the Lorentz factor above 10, the 
mechanism still can work efficiently, but only at favourable con- 
ditions: a soft spectrum of external isotropic radiation and a weak 
(dynamically not important) magnetic field. When the Lorentz fac- 
tor is below ~ 8 or the AGN luminosity is below 10 43 i?i7erg s _1 
the photon breeding mechanism does not work. 

The mechanism can work in different astrophysical contexts: 
in the immediate vicinity of the accretion disc, in the broad line re- 
gion (probably, most efficiently), and in the IR radiation field pro- 
duced by the dust at a parsec scale. There could be other sites for 
operation of the photon breeding mechanism (e.g. stellar radiation 
field at a kpc scale, cosmic microwave background radiation at 100- 
kpc scale). 

There exist a clear signature of the mechanism which can be 
observed by GLAST in the spectra of moderately bright blazars: a 
synchrotron cutoff at a few GeV range. 

There exist some problem with the straightforward interpre- 
tation of the blazar spectra. The low-energy synchrotron peak re- 
produced by our simulations is less prominent than the one ob- 
served in the spectra of gamma-ray quasars. An additional mecha- 
nism of gradual pair reheating is necessary. This could be associ- 
ated with a diffusive Fermi acceleration of pairs, produced by the 
photon breeding, by fluid perturbations. Such reheating mechanism 
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should have a moderate energy budget and seems very natural af- 
ter a strong impact of radiation on the jet. Importantly, the photon 
breeding solves the problem of injection of relativistic electrons for 
Fermi acceleration: they are already produced by high-energy pho- 
tons. 

Fast deceleration of the outer jet layers results in the fast 
spine-slower sheath structure which is also implied by various ob- 
servations. Particularly, the differential deceleration can resolve the 
contradiction between a requirement of the large Doppler factor in 
the gamma-ray emitting region and the low Lorentz factor of jets in 
BL Lacs derived from the radio and optical observations. 

Despite the fact that the photon breeding mechanism has not 
been known before and may seem somewhat exotic, it is actually 
much simpler in the description than the diffusive Fermi accelera- 
tion, which depends on complicated plasma phenomena. Once the 
velocity pattern for the fluid, the magnetic field and the external ra- 
diation field are specified, the fate of each high-energy photon and 
its descendants can be reproduced (in statistical sense) from first 
principles, because the interaction cross-section are known with 
high accuracy. The question whether a photon produces a runaway 
avalanche or not, can be answered exactly. However, there still re- 
main some details which require a more elaborated treatment of the 
fluid dynamics and a more complete model for the jet. Particularly, 
it would be important to study the jet launch, the photon breeding, 
and the jet deceleration with the resulting gamma-ray emission in 
one numerical experiment. 
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